ON A CURIOUS PROPERTY OF BELL NUMBERS 



Zhi-Wei Sun and Don Zagier 

Abstract. In this paper we derive congruences expressing Bell numbers 
and derangement numbers in terms of each other modulo any prime. 



1. Introduction 

Let B n denote the nth Bell number, defined as the number of partitions 
of a set of cardinality n (with B = 1). In 1933 Touchard [T] proved that 
for any prime p we have 

B p+n = B n + B n+ i (mod p) for all n = 0, 1, 2, (1) 

Thus it is natural to look at the numbers B n (mod p) for n < p. In [S] , the 
first author discovered experimentally that for a fixed positive integer m 
the sum Yln=o (~ m ) n modulo a prime p not dividing m is independent 
of the prime p, a typical case being 

p-i B 

7 — <J^T ^ —1853 (mod p) for all primes p^2. 

n=0 ^ ' 

In this note we will prove this fact and give some related results. 

Our theorem involves another combinatorial quantity, the derangement 
number D n , defined either as the number of fixed-point-free permutations 
of a set of cardinality n (with Dq = 1) or by the explicit formula 

^ = t'^T (" = 0,1,2,. ..)• (2) 

k=0 
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Theorem 1. For every positive integer m and any prime p not dividing 
m we have 

E 7^ = ("ir-^m-i (mod p) . (3) 



Using Eo<m< P (" m )" k = -$n,k (mod p) for fc, n G {1, . . . , p - 1}, we 
immediately obtain a dual formula for B n (n < p) in terms of Do, . . . , D v -i- 

Corollary. Let p be any prime. Then for all n = 1, . . . ,p — 1 we have 

p-i 

B n =J2 (-l) m D m -i (-m) n (mod p) . 

m=l 

For the reader's convenience we give a small table of values of B n and D n . 



n 





1 


2 


3 


4 


5 


6 


7 


8 




1 


1 


2 


5 


15 


52 


203 


877 


4140 




1 





1 


2 


9 


44 


265 


1854 


14833 



We will show Theorem 1 in the next section and derive an extension of 
Theorem 1 in Section 3. 

2. Proof of Theorem 1 

We first observe that it suffices to prove (3) for < m < p, since both 
sides are periodic in m (mod p) with period p. For the left-hand side this 
is obvious and for the right-hand side it follows from (2), which gives the 
expression (—l) n D n = J2^ (— l) r n(n — 1) ■ • • (n — r + 1) for D n (mod p) 
as a terminating infinite series of polynomials in n. 

We will prove (3) for < m < p by induction on m. Denote by S m the 
sum on the left-hand side of (3) , where we consider the prime p as fixed and 
omit it from the notation. Since D n = nD„_i + (— l) n for n = 1, 2, 3, . . . 
(obvious from (2)), we have to prove the two formulas 

Si = 1 (mod p) , mS m = Si - S m+ i (mod p) . (4) 
Recall that the Bell numbers can be given by the generating function 
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equivalent to the well-known closed formula 



oo 

e — ' r! 

r=0 



Since the function y = e e * 1 satisfies y' = e x y, this also gives the recursive 
definition 

S = 1 , B n+1 = J2 (l) B k for all n > . (6) 

k=o ^ ' 

This recursion is the key ingredient in proving (4). 

For the first formula in (4) we use (6) with n = p — 1 to obtain 

Si = y7-l) fc B fc = V f. B fe = B p -B (modp), 
fc=i fc=i V k J 

so it suffices to prove that B p = 2 (mod p). This is a special case of 
Touchard's congruence (1), but can also be seen by writing (5) in the form 

n=0 Kr<p ^ 

to get Bp/p! = l/p\ + (p-integral) + l/p\. 
Now using Fermat's little theorem we have 

-mS m = E(- m ) P " 1_n EQ^ 

n=0 fc=0 ^ ' 

p-2 p-k-2 , , s 

ee J>l) fe B fe £ P J V*- 1 -' (r = n-k) 

fc=0 r=0 ^ ' 

p-2 

= J2(- 1 ) kB k((m+l) p - 1 - k -l) = S m+1 - S 1 (modp) 

fc=0 

for 1 < m < p — 2. This completes the proof of (4) and the theorem. 

3. An extension of Theorem 1 

Recall that for nonnegative integers n and k the Stirling number S(n, k) 
of the second kind is the number of ways to partition a set of n elements 
into k groups. Obviously 

n 

B n = ^S(n,k). 

fc=0 
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The Touchard polynomial T n (x) of degree n is given by 

n 

T n (x) = ^S(n,k)x k . (7) 



k=0 



Note that T n (l) = B n . Similar to the recursion for Bell numbers, we 
have the recursion 



T n+1 (x) = xJ2(f)Tk(x). 

k=0 ^ ' 



(8) 



Let p be a prime and let Z p be the ring of p-adic integers. For two 
polynomials P(x),Q(x) G 1> p [x], by P(x) = Q(x) (mod p) we mean that 
the corresponding coefficients of P(x) and Q(x) are congruent modulo p. 

Our next theorem is a further generalization of Theorem 1. 

Theorem 2. For every positive integer m, we have 

("*>" E ^ - -^£<^(-*)< (modp) (9) 

0<n<p ^ ' 1=0 

for any prime p not dividing m. 

As a consequence, if x is a p-adic integer not divisible by p, then 

E^ -m n -E) m E^ fc v ; v 1 ' y ' 

0<n<p y J v ' k=0 

In particular, 



E , n . = (mod p) for p ^ 2, 

(-2) n x v y 



0<n<p 



^ T n Qr) - 2x + 2 

2^ TE^ET = ^ (modp) forp^3, 



0<n<p 

T n (x) — 3x 2 + 6x — 6 



(Ei^T = ~3 ( mod ^) forp^2. 



0<n<p 



Although we can show Theorem 2 via a slight modification of the proof 
of Theorem 1, below we prove Theorem 2 by a new approach. 
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Proof of Theorem 2. Observe that 

Tnjx) = P ^ Ek=iS(n,k)x k = ^ k P ^ S(n,k) 

It is known that 

S ( n > ^ = ^ E (*) (- 1 ) fc "^' n for a11 M = 0, 1, 2, . . . . 

' 3=0 ^ 7 

Thus 

1 (-m) n f-' kl ^-rf-f V ?/ V m 



fe=l n=l j=l 

For each j e {1, . . . ,p — 1}, if p | m + j then 
p-i / • \ n p-i 



s E 1 s - 1 (modJ,) - 



n=l v n=l 

if p \ m + j then 

p-l / . \ n P-I 



n=l v 7 n=0 v 7 J/ 



with the help of Fermat's little theorem. 

Let r denote the least positive residue of — m modulo p. By the above, 

(_ TO )n Z^yAU 1 J 2^1 m 

n=l v 7 fc=l 3 = 1 XJ 7 n=l v 

^Ei(!)(-r(-.) = t^EgM ? ). 

fc=l V 7 fc=r V ' 



Therefore 



(_„)-. T = F ( ~ J> , modp) 

v ; ^ (-m) n r! ^ (p + Z-m-r)! 

n=l v 7 l=m+i — p v 7 

So it remains to show that 

(m — 1)! f (— l) r+1 /(r!(/j + / — m — r)!) (mod p) if m + r — p ^ / < m, 



11 [0 (mod p) if ^ / < m + r — p. 
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IfO^Z<m + r — p, then we have (m — = (mod p) since 

I < m + r — p < m and m + r = (mod p) . 
Now suppose that m + r — p ^ I < m. Then 

(-l) r+1 fp + l-m 



(p + I — m)\ \ r 

( _ 1)r +i uTJ^iP-s) fi 



m 

(p — 1)1 \ r 

_ f^m-i-ifa _i_iy( m - l + r ~ l \ ( by Wilson's theorem) 

\ m — l — 1 J 

= (-l)™-'-^™- /-!)!( J = (m ~ 1)! (modp) (asp | m + r) 



m-l-1) l\ 

In view of the above, we have completed the proof of Theorem 2. □ 
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